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1  L  _ _ 

ON  THE  POSSIBILITY  OF  HIGH  P0V\T:R  GYROTRONS 
FOR  SUPER  RANGE  RESOLUTION  RADAR 
AND  ATMOSPHERIC  SENSING 


1,  Introduction  and  Review  of  NRL  Cryotrons  at  About  100  GHz 

High  power  gyrotrons  have  been  develojjed  over  the  last  15 
years  in  the  United  Statesl'5,  the  Soviet  Union^.^,  Europe8.9,  Japan^O, 
and  China^^.  The  research  has  principally  been  directed  towards 
developing  sources  for  electron  cyclotron  heating  of  fusion  plasmas. 
The  emphasis  is  on  developing  cw  sources  with  high  power. 
Commercial  gyrotrons  are  available  from  Varian  Associates,  Palo 
Alto,  Cal,  with  300  kW  at  28  GHz,  200  kW  at  60  GHz,  and  100  kW  at 
140  GHz.  Varian  is  now  working  with  a  whispering  gallery  mode 
gyrotron,  originally  demonstrated  at  MIT,  to  develop  a  1  MW 
gyrotron  at  140  GHz^2  the  Naval  Research  Laboratory,  a  94  GHz 
gyrotron  has  been  developed  at  150  kWl3.l4^  anj  the  quasi-optical 
gyrotron  has  demonstrated  600  kW  of  peak  power  and  tunability 
from  about  80-130  GHz ^ 5-21  This  power  and  tunability  is  a  very 
new  capability,  and  we  are  looking  for  additional  applications.  We 
examine  here  several  possible  applications  involving  high  range 
resolution  radar  and  atmospheric  sensing.  Since  the  sources  we 
discuss  are  oscillators  rather  than  amplifiers,  all  of  the  applications 
depend  on  a  variant  of  coherent  on  receive  data  processing  in  radar 
systems.  Specifically  the  entire  amplitude  and  phase  modulation  of 
the  transmitted  and  received  waveforms  are  recorded  digitally  and 
then,  processed,  most  likely  off  line.  The  ability  to  do  this  results 
from  the  extremely  rapid  advances  in  data  processing  technology. 

For  instance  a  500  MHz-1  Ghz  A  to  D  converter,  which  is  now  state  of 
the  art,  did  not  exist  a  few  years  ago.  Similarly,  if  data  processing 
capabilities  increase  at  the  same  rate  in  the  future,  a  few  years  from 
now  real  time  data  processing  could  be  utilized  for  these  applications. 
Section  2  will  discuss  in  some  detail  the  data  processing  techniques. 

A.  The  NRL  94  GHz  Gyrotron 

The  94  GHz  gyrotron  is  a  TE13  mode  gyrotron  which  has 
generated  150  kW  of  power  in  a  Ips  pulse  at  an  efficiency  of  15- 
20%.  With  appropriate  cooling  on  the  collector  and  window,  it  could 
operate  at  an  average  power  of  10  kW.  This  is  in  fact  only  10%  of 
the  average  power  Varian  has  demonstrated  in  their  140  GHz 
gyrotron.  A  mode  convertor  has  been  developed  by  General  Atomics, 
San  Diego,  Cal,  and  is  currently  in  place  on  our  gyrotron.  Conversion 
to  the  TE]  1  mode  has  been  demonstrated.  Figure  1  is  a  photograph  of 
the  gyrotron  laboratory.  The  equipment  is  compact  and  can  easily  be 
transported  on  a  small  truck.  This  capability  to  radiate  average 
powers  of  order  10  KW  at  94  GHz  inu,  the  atmosphere  represents  a 
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very  new  capability.  By  contrast,  the  highest  source  available  now  is 
the  Lincoln  Laboratory  Kwajalein  radar  tube  which  has  a  peak  power 
of  6  KW  and  an  average  power  of  about  600  Watts  and  a  bandwidth 
of  about  1  GHz22. 


B.  The  NRL  Quasi -Optical  Gyrotron 

The  quasi-optical  gyrotron  is  a  very  different  sort  cf  a  gyrotron 
in  an  optical  resonator.  It  has  been  developed  at  both  NRL  and 
Lausanne23.24  beam  propagation  and  the  magnetic  field  are 

both  transverse  to  the  axis  of  the  resonator.  The  mode  output  is 
taken  by  diffraction  around  the  mirror.  The  diffraction  is 
characterized  by  a  fractional  loss  v  per  round  trip  of  the  radiation. 
Typically  our  experiment  runs  with  v's  between  1  and  5%. 

Competition  between  different  longitudinal  modes  is  an  important 
factor  in  the  operation  of  the  device.  If  the  mirrors  are  very  close 
together  the  mode  density  is  low  and  it  oscillates  in  oniy  one  mode  at 
a  time.  However  if  the  mirrors  are  further  apart,  it  typically 
oscillates  in  several  distinct  longitudinal  modes.  Figure  2  is  a 
photograph  of  the  quasi-optical  gyrotron.  It  is  about  as  portable  as 
the  conventional  gyrotron. 

At  NRL,  two  separate  sets  of  experiments  on  the  configuration 
with  wide  mirror  separation  have  been  conducted;  one  used  a  Varian 
VUW  8010  electron  gun  that  can  run  up  to  about  75  kV  and  20 
Amps,  and  another  used  a  Varian  VUW  8144  electron  gun  which  can 
operate  up  to  about  100  kV  and  60  Amps.  Figure  3  is  a  plot  of 
output  power  and  efficiency  for  the  latter  gun,  as  a  function  of 
current,  at  a  center  frequency  of  about  120  GHz.  The  power  peaks  at 
600  kW  and  the  efficiency  peaks  at  12%. 

One  of  the  most  attractive  features  of  the  quasi-optical 
gyrotron  is  its  tunability.  As  the  frequency  changes,  the  only  thing 
that  changes  in  the  cavity  is  the  output  coupling,  which  slowly 
decreases  as  a  function  of  frequency.  The  quasi-optical  gyrotron  can 
be  tuned  in  three  ways;  •  The  magnetic  field  can  be  varied  over  a 
fairly  wide  range.  The  power  is  relatively  insensitive  to  field.  Figure 

4  shows  a  plot  of  mode  frequency  and  power  as  a  function  of 
magnetic  field  for  the  lower  power  electron  gun.  The  device  can  tune 
over  about  half  an  octave  from  about  80  to  130  GHz.  •  The  device 
can  be  tuned  with  the  beam  voltage.  Since  the  frequency  scales  as 
B/y,  the  frequency  can  be  raised  by  by  lowering  the  voltage.  Figure 

5  shows  the  Voltage  tuning  curve  for  the  low  power  gun  experiment. 


input  power  is  proportional  to  it.  •  The  cavity  can  be  tuned  by 
adjusting  the  mirror  separation,  which  can  be  done  without  breaking 
vacuum. 

It  is  also  possible  that  the  quasi-optical  gyrotron  could  operate 
at  the  second  harmonic,  although  at  lower  power.  An  experiment  can 
be  designed  with  the  harmonic  having  a  lower  starting  current  by 
exploiting  the  fact  that  the  radiation  beam  has  lower  output 
coupling25.26.  The  Lausanne  experiment  has  already  seen  evidence 
of  harmonic  operation^^.  If  the  quasi-optical  gyrotron  can  run  with 
harmonic  operation,  it  opens  the  possibility  of  a  single  source 
covering  the  entire  spectral  region  from  80  to  260  GHz.  This  would 
provide  about  1  kW  of  average  power  in  fundamental  mode 
operation  and  about  10-100  Watts  of  average  power  in  harmonic 
operation. 


C.  Applications 

With  this  review  of  the  characteristics  of  the  NRL  gyrotrons,  we 
continue  with  a  discussion  of  potential  applications  for  fielded 
systems.  In  any  outdoor  millimeter  wave  system,  the  question  of 
interaction  of  the  millimeter  waves  with  the  atmosphere  is  a  crucial 
question.  The  propagation  of  millimeter  wave  power  in  the  air  has 
been  studied  extensively  by  Liebe,27  and  we  use  his  model.  Figure  6, 
provided  by  Liebe,  shows  the  results  of  his  model  for  attenuation 
and  phase  shift  at  sea  level  at  25°C  as  a  function  of  frequency  and 
relative  humidity.  Also  a  very  useful  information  source  is  the  Near 
Millimeter  Wave  Technology  Base  Study,  edited  by  Kulpa  and 
Brownes. 

Possible  applications  for  wideband,  high-power  millimeter 
waves  include  very  wide  band  radar  operation  and  atmospheric 
sensing.  There  is  now  interest  in  ultra  wideband  radars  for  fielded 
systems29.30  These  ultra  wideband  systems  are  proposed  with 
bandwidths  from  1-3  GHz.  However  the  quasi-optical  gyrotron,  in  its 
fundamental  mode  operation  has  already  demonstrated  more  than 
40  GHz  of  tunability.  Section  2  shows  how  this  tunability  can 
translate  to  bandwidth,  giving  2-3  mm  range  resolution. 

Section  3  discusses  the  possibility  of  using  millimeter  wave 
gyrotrons  to  obtain  radar  images  of  clouds.  Since  clouds  are  made 
from  aerosols  whose  diameter  is  much  less  than  a  wavelength,  the 
scattering  is  in  the  Rayleigh  regime.  The  cross  section  scales  as  a 
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scattering  is  in  the  Rayleigh  regime.  The  cross  section  scales  as 
giving  tremendous  potential  advantage  to  short  wavelength  systems. 
These  systems  are  restricted  principally  by  the  unavailability  of  high 
power  sources.  In  fact  a  series  of  radar  images  of  clouds  has  been 
obtained  by  Lhermitte3l-34  by  using  an  extended  interaction 
oscillator  at  94  GHz  with  an  average  power  of  4  W.  We  discuss  how 
the  use  of  rf  sources  with  much  higher  average  power  and  tunability 
can  greatly  extend  this  capability  by  using  standard  meteorological 
data  processing  schemes35,36  but  applied  to  oscillators. 

Section  4  discusses  the  use  of  high-power  millimeter-wave 
sources  to  remotely  sense  the  characteristics  of  clear  air  turbulent 
flow.  At  longer  wavelengths,  radar  backscatter  from  clear  air 
turbulence  is  routinely  used  to  measure  average  wind  velocity37-40. 
However  very  high  power  and  very  large  antenna  size  are  necessary 
and  these  measurements  are  almost  always  done  from  fixed  cites. 

(Or  else  the  experiments  are  done  at  very  short  range.)  For  instance 
the  antenna  in  Ref.  39  was  26  m  in  diameter  and  the  wavelength  was 
about  a  foot.  The  peak  power  was  about  4  MW.  Thus  the  scattering 
was  from  fairly  long  wavelength  components  of  the  turbulent 
spectrum.  The  availability  of  high  power  millimeter  wave  sources 
opens  up  the  possibility  of  scattering  from  much  shorter  wavelength 
fluctuations  with  a  moveable  facility.  Specifically  we  investigate  the 
way  the  minimum  scale  length  of  the  turbulent  flow,  which  is 
directly  correlated  with  the  turbulent  dissipation,  can  be  detected. 

Section  5  discusses  the  remote  measurement  of  humidity  by 
two  techniques.  The  first  is  a  line  of  sight  propagation 
measurement'll,  but  now  over  a  long  horizontal  path.  The  second 
measures  humidity  profile  in  clouds  by  comparing  radar  scatter  at 
two  frequencies.  Section  6  examines  the  analogous  application  of  the 
tunability  of  the  quasi-optical  gyrotron  to  the  measurement  of  trace 
impurities  in  the  upper  atmosphere  by  two  way  propagation 
measurements  (ie  reflection  from  an  orbiting  satellite). 
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2.  The  Quasi-Optical  Gyrotron  for  a  High-Power,  Ultra  Wide  Band  Radar 

In  this  section  we  show  how  an  oscillator  based  radar  can  achieve  a  range 
resolution  equal  to  c/2Af.  In  the  quasi-optical  gyrotron,  by  stepping 
through  a  sequence  of  modes,  Af  can  be  made  equal  to  40  GHz,  the  full 
tuning  range.  The  key  to  oscillator  use  is  in  the  data  processing  in  the 
reciever.  The  transmitted  and  received  pulse  are  each  digitally  recorded 
and  an  effective  matched  filter  is  generated  numerically. 

A.  Radar  Range  Equation 

The  power  and  wavelength  of  current  gyrotrons  make  the  radar 
resolution  of  very  small  targets  possible  at  very  large  range.  The  radar 
equation  for  the  range  R  is 

R4  exp  2aR  =  A,2G2Pta/(47c)3Pr. 

where  a  is  the  atmospheric  attenuation,  X  is  the  wavelength,  G  is  the 
antenna  gain,  Pt  is  the  transmitted  power,  o  is  the  target  cross  section,  and 
Pr  is  the  received  power.  The  antenna  gain  is  related  to  wavelength  and 
diameter  d  of  an  ideal  radiator  by  G  =  We  assume  that  for 

reasonable  radar  performance  the  received  power  is  100  times  the  noise 
power,  Pn  =  kTAf,  where  Af  is  the  receiver  bandwidth  and  the  temperature 
is  assumed  to  be  SOO^K.  We  consider  two  possibiUties  for  the  bandwidth. 
We  first  assume  that  the  bandwidth  is  the  width  of  the  individual  line  of  the 
quasi-optical  gyrotron,  (o/Q,  about  15  meters  if  Q  =  ICH.  In  the  second 
case,  we  consider  the  bandwidth  to  be  pulse-limited;  that  is  for  a  10  ps 
pulse  the  bandwidth  is  10^,  so  range  resolution  in  an  individual  pulse  is  15 
km.  There  is  a  wide  frequency  window,  from  about  70-1 10  GHz  where 
the  atmospheric  attenuation  in  clear  weatlier  is  less  than  0.5  dB/km  at 
relative  humidity  below  50%.  Then  if  a  one  meter  dish  is  used  as  the 
antenna  and  the  transmit  power  is  300  kW,  the  maximum  range  of  a  0.1 
cm2  target  is  15  km  for  the  wide  band  case  and  is  30  km  for  the  narrow 
band  case. 


B.  Single  Pulse  Radar  Return 

For  a  target  and  transmitter  that  are  both  stationary,  there  is  no 
frequency  shift  upon  reflection.  Therefore  since  the  transmitted  signal  is 
square  integrable  in  time  from  minus  infinity  to  plus  infinity,  it  can  be 
written  as  a  Fourier  transform. 
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Et(t)  =  J  df  Gt(f)  exp  27cift  =  et(t)  cos[27tfOt  +  (t)t(t)].  (1) 

Here  fo  is  assumed  to  be  the  carrier  frequency  and  e  and  (j)  are  assumed  to 
vary  very  slowly  in  time  compared  with  fo.  We  adopt  a  notation  using  E 
as  a  dependent  variable  in  the  time  domain  and  G  as  a  dependent  variable  in 
the  frequency  domain.  We  assume  then  that  e  and  (j)  can  be  digitized  and 
computer  analyzed.  The  digitizing  does  not  have  to  be  nearly  fast  enough 
to  resolve  fo.  We  discuss  later  just  what  the  requirements  on  the  digitizer 
must  be. 

The  fundamental  assumption  that  we  make  here  is  that  a  return  signal 
can  be  associated  with  one  and  only  one  transmitted  pulse.  Particularly, 
there  does  not  appear  to  be  any  way  that  a  radar  based  on  an  oscillator 
signal  can  separate  returns  from  different  range  regions.  Thus  for  a  pulse 
separation  Xp,  returns  from  signals  at  distance  greater  than  cXp/2  cannot  be 
processed  correctly.  We  now  consider  the  returned  signal  associated  with  a 
particular  transmitted  signal. 

The  returned  signal  can  be  written  as 

Er(t)  =  J  df  Gr(f)  exp  27cift  s  J  df  s(f)  Gt(f )  exp  27rift  (2) 

H  er(t)  cos  [27cfot  +  -  IiSiet(ti)cos[27tfoti  +  (l)t(ti)  ], 

where  s  is  the  coefficient  relating  the  transmitted  to  returned  signal  at 
frequency  f.  Here  s  =  ^GVaexp(-aR)/(47c)3/2R2  jn  the  last  line  of 
Eq.(2),  we  assume  that  the  returned  signal  is  a  summation  over  discrete 
scatterers  at  r2mge  Ri.  The  quantity  U  is  t-2Ri/c. 

The  key  to  using  an  oscillator  rather  than  an  amplifier  then  is  to 
digitally  record  Et(t),  or  equivalently  et(t)  and  (|)t(t),  on  each  pulse  and 
calculate  Gt(f)  by  taking  a  Fourier  transform.  That  is  a  very  fast  operation 
on  a  computer  using  fast  Fourrier  transform  routines.  Then  the  same  thing 
can  be  done  on  the  returned  signal  to  get  Gr(f).  This  then  allows  s  to  be 
calculated  as  a  function  of  frequency  as 

s(f)  =  Gr(f)/Gt(f).  (3) 

In  practice,  the  transmitted  signal  does  not  have  a  meaningful  reciprocal 
outside  its  bandwidth.  Thus  the  usable  bandwidth  for  an  oscillator  based 
system  is  less  than  the  total  bandwidth  of  the  transmitted  signal. 
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From  the  assumption  that  the  system  is  stationary  in  time,  we  can 
calculate  the  radar  return  had  the  transmitted  pulse  been  a  desired  pulse. 
Using  a  subscript  d  to  denote  the  desired  pulse  on  either  the  transmitted  or 
returned  pulse,  we  find 


Grd(f)  =  s(f)  Gtd(f).  (4) 

The  desired  pulse  has  knoivn  properties,  which  we  specify.  Typically  the 
desired  pulse  is  a  Gaussian  with  linear  frequency  chirp. 

To  get  e(t)  and  (t)(t),  we  send  the  signal  and  two  reference  signals, 
proportional  to  sin  and  cos  of  27cfot  through  a  mixer  followed  by  a  low 
pass  filter.  This  will  give  two  signals  proportional  to  e(t)sin(l)(t)  and 
e(t)cos(()(t),  from  which  e  and  <()  can  be  obtained.  They  both  vary  very 
slowly  in  time  compared  with  fo  and  can  be  sampled  and  digitized.  For 
instance,  a  500  MHz  A/D  convertor  with  8  bit  digitizer  is  now  state  of  the 
art.  It  will  let  the  signal  be  digitized  at  this  rate  and  allow  256  values  of 
the  function  on  a  linear  or  logarithmic  scale. 

Of  course  the  digitized  samples  must  be  processed  by  a  discrete 
Fourier  transform  (DFT).  To  illustrate  the  nature  of  the  range  resolution 
obtained,  the  discrete  calculations  are  approximated  by  continuous  signals. 
Now  let  us  Fourier  transform  the  continuous  transmitted  signal.  It  is 

Gt(f)  =  /  dt  et(t)  cos[27tfot  +  <|)t(t)]  exp(-2Kift).  (5) 

Assuming  that  the  fastest  variation  is  from  fQ  and  that  we  are  looking  at  f  > 
0  (the  assumption  is  always  that  fQ  >  0),  we  find  that 

Gt(f)  =  0.5j  dt  et(t)  exp  12Ki5ft+i<J)i(t)],  f  >  0  (6) 

where  6f  =  fo-f.  Thus  the  Fourier  transform  can  be  calculated  from  only 
a  knowledge  of  e  and  ')  and  only  by  operating  on  a  time  scale  much  slower 
than  fo,  that  is  a  time  scale  on  which  the  data  can  be  digitized.  From  a 
calculation  of  both  Gt(f)  and  Gr(0»  the  cross  section  s(f)  can  be  calculated. 
From  this  we  can  calculate  Grd(f)  from  the  knowledge  Gtd(f). 

The  next  step  is  calculating  of  erd(t)  and  (t)rd(t).  From  the  definition 
of  the  Fourier  transform,  we  find  with  a  simple  calculation  that 


oo 

erd(t)  exp  i(t)rd(t)  =  lodf  Grd(0  exp  (-27ii5ft).  (7) 
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Thus  Crd  (t)  and  (})rd(t)  can  be  obtained  from  the  returned  signal  entirely 
by  considering  the  slow  time  scale  and  without  digitizing  anything  at 
frequency  fQ. 

Once  we  have  the  e's  and  (j)’s  for  the  desired  transmitted  and 
returned  signals,  we  form  the  matched  filter  response  for  the  n^^  range  bin, 

Mn  =  J  dt  etd(t-tn)erd(t)  exp  i[(t)rd(t)  -  <l)td(t-tn)]  = 

XiJdt  Si  etd(t-tn)  etd(ti)  exp  i[(l)td(ti)  -  (l)td(t-tn)  -47ifoRi/c].  (8) 

Note  that  with  this  sort  of  signal  processing,  there  is  no  single  matched 
filter  for  the  radar;  rather,  a  separate  matched  filter  is  generated 
numerically  for  each  radar  pulse.  This  filter  response  then  characterizes 
the  cross  section  of  the  range  interval  centered  at  t=tn  or  Rn  =  ctn/2.  If  we 
assume  that  the  desired  transmitted  signal  is  a  Gaussian  envelope,  etd(t)  = 
eoexp-(t/T)2,  with  a  linear  frequency  chirp,  (})td  =  O.Sat^,  the  integrals  can 
be  done  for  both  the  matched  filter  response  and  for  the  Fourier 
transforms.  The  results  are 

Gt(f)  =  {^7teoT/[l-iaT2/2]I/2}  exp['7i25f2(x-2+i(xy^X-'^  +  a2/4j]  (9) 
and 

Mn  =  (V7c/2)eo^T  Xi  Siexp(-47tifoRi/c) 

X  exp{-(tn-2Ri/c)2[l/2T2  +  a2T2/8]}.  (10) 

The  range  resolution  is  determined  by  the  fall  off  of  the  Gaussian.  If  this 
is  dominated  by  the  frequency  chirp,  the  range  resolution  goes  roughly  as 
caT/V8  =  c/Af  where  Af  is  the  frequency  spread  of  the  transmitted  pulse. 


C.  Data  Processing  for  Time  Limited/  Band  Limited  Functions 

Although  a  ftinction  that  is  limited  in  both  time  and  band  is  not 
possible,  it  is  a  useful  approximation  and  idealization  and  gives  important 
insight  into  the  needs  of  the  data  processing  system.  The  Gaussian 
waveform  with  the  linear  chirp  is  a  reasonable  approximation  to  a  time- 
limited,  band-limited  function,  and  other  functions  may  be  better  still.  The 
transmitted  pulse  is  assumed  to  be  limited  to  time  between  0  and  Tt.  The 
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received  pulse  is  then  limited  to  times  between  2Rmin/c  and  Tt+2Rmax/c, 
where  Rmin  and  R  max  are  the  minimum  and  maximum  range.  Depending 
on  the  spatial  extent  of  the  target,  the  time  of  the  returned  pulse  can  be 
comparable  to  or  much  longer  than  the  transmitted  pulse.  Note  that  if  the 
returned  signal  is  much  longer  than  the  transmitted  signal,  more  spectral 
density  is  needed  to  resolve  it  in  the  Fourier  domain.  However  the 
maximum  bandwidth  required  is  determined  by  the  transmitted  signal. 

Since  the  incident  and  returned  signals  are  limited  in  time  to  a  time 
we  can  To,  we  can  regard  these  signals  as  periodic  with  period  To  as  long 
as  we  are  careful  to  always  stay  within  the  time  0  <  t  <To,  Thus  the  signal 
can  be  regarded  as  a  discrete  Fourier  summation 

f(t)  =  Xn  g(n)  exp  27cint/To,  0  <  t  <  To.  (11) 

If  the  signal  is  band-limited  as  well  as  time-limited,  then  g(n)  is  non  zero 
only  for  Inl  <  N.  Then  f(t)  is  needed  only  at  2N  discrete,  equally  spaced 
values  of  t  to  solve  for  the  2N  values  of  f(n).  Thus  the  waveforms  must  be 
sampled  twice  per  period  at  the  maximum  frequency,  and  this  gives  all  the 
information  there  is  to  be  had  about  the  function. 


D.  H’qh  Range  Resolution 

Now  let  us  consider  using  the  full  bandwidth  of  the  NRL  quasi- 
optical  gyrotron  in  a  stepped  frequency  mode  to  achieve  very  high  range 
resolution  of  ar.  isolated,  stationary  target.  Let  \i^  consider  the  target  to  be 
near  the  center  of  a  range  cell  for  a  single  pulse.  Then,  according  to 
Eq.(ll), 


M(q)  =  (V7c/2)eo“T  Xi  Si  exp  (-47uifqRi/c),  ( 1 2) 

where  now  we  use  an  index  q  for  the  mode  frequency  (fq=  qc/2L=  qw)  to 
denote  that  it  will  change  from  pulse  to  pulse  and  we  have  deleted  the  index 
n  on  M  because  everything  is  assumed  to  take  place  in  a  single  range  cell, 
that  is,  at  a  single  value  of  n.  Here  fq  is  the  frequency  of  the  local 
oscillator,  which  is  nearly  equal  to  the  frequency  of  the  quasi-optical 
gyrotron  when  it  operates  in  ♦he  desired  mode.  Also  we  will  use  a 
dimensionless  range  to  the  target  pi  =  2wRi/c.  As  the  quasi-optical  gyrotron 
shifts  from  one  mode  to  another,  q  will  increase  from  a  minimum  value  qi 
up  to  a  maximum  value  q2- 
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Note  that  the  reference  frequency  fq  must  be  detemined  very 
accurately  because  any  error  in  4jcfqR/c  must  be  significantly  less  than  one 
from  pulse  to  pulse.  Thus  as  one  goes  from  frequency  to  frequency,  the 
local  oscillator  frequency  must  be  known  very  accurately  (but  not 
necessarily  the  high-power  oscUlator  frequency).  Fortunately  such 
sources  have  been  developed.  For  instance  the  Hughes  Series  4774xH  phase 
locked  Gunn  oscillator  seems  to  have  sufficient  stability.  The  phase  is 
locked  to  a  high  harmonic  of  a  low  frequency  source  (tens  of  megahertz) 
and  the  Hughes  catalog  shows  that  phase  noise  is  down  60db  compared  to 
the  main  signal  10  Hz  away  from  Ae  carrier  frequency,  and  the  amplitude 
noise  is  even  lower. 

For  now,  assume  that  all  ranges  pi  are  known  and  form  the 
summation 

q2  q2 

Zql  M(q)  exp  27cipjq  =  (V7c/2)eo2T  SiZql  Si  exp  -27ciq(pi-pj) 

=Si  Si  [exp  -7ti(pi-pj)(q2-qi)]{simi(pi-pj)(q2+l-oi)/sin7t(pi-pj)}.  (13) 

The  last  line  in  Eq.(14)  comes  from  summing  up  the  geometric  series  in  q. 
The  term  in  the  brackets  is  very  much  like  a  delta  function  and  maximizes 
at  pi=pj.  Thus,  regarding  the  left  hand  side  of  Eq,(14),  which  is  a  function 
of  pj,  we  see  that  contributions  to  it  come  only  from  nearby  targets.  The 
contribution  to  the  left  hand  side  of  Eq.(14)  then  comes  from  targets  within 
the  first  maximum  of  the  sin  function  in  the  brackets,  or  Rj-Ri  <  c/2(f2-fi), 
where  we  have  gone  back  to  dimensional  units.  This  is  the  conventional 
for  range  resolution  a  radar  of  having  frequency  spread  f2-fl. 

Let  us  now  approximate  the  target  as  a  series  of  discrete  scatterers  at 
ranges  given  by  Pp  =  p/(q2+l-qi)»  where  p  is  an  integer  that  ranges  from  a 
minimum  value  of  pi  to  a  maximum  value  of  p2.  Then  the  brackets  on  the 
right  hand  side  of  ^.(14)  becomes  proportional  to  a  Kronecker  delta  and 
(changing  the  indices  of  summation  from  ij  to  pp’)  we  find 

q2 

Iql  M(,q)  exp  27tipq/(q2+l -qi)  =  (q2+l -qi)(P2+l -pi)Ctp.  (14) 

Thus  if  the  targets  are  assumed  to  be  discrete  points,  separated  by  the  range 
resolution  and  localized  at  a  grid  in  range  that  is  reciprocal  to  the  grid  of 
discrete  frequencies,  the  cross  section  at  each  position  is  simply  the  discrete 
Fourier  transform  of  the  filter  response  as  a  function  of  q.  The  maximum 
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number  of  elements  that  can  be  resolved  is  then  q2  -  qi  +1  as  long  as  p2-pi 

=  q2  -qi. 

Notice  however  that  if  p  increases  by  q2+l-ql.  the  left  hand  side  of 
Eq.(14)  is  unchanged.  Thus  there  is  an  ambiguity  if  the  spatial  extent  of 
the  target  is  larger  than  L,  which  is  the  mirror  separation  in  the  quasi- 
optical  gyrotron.  This  ambiguity  can  be  mitigated  by  varying  the  mirror 
separation. 

Figure  7  is  a  block  diagram  of  the  circuitry,  for  the  case  of 
operation  in  one  mode  at  a  time.  Notice  that  after  the  A/D  converter,  the 
blocks  do  not  correspond  to  physical  circuit  components,  but  to  elements  or 
subroutines  of  a  computer  code. 

E.  Rapid  Range  Imaging  of  Isolated  Targets 

We  have  seen  how  the  quasi-optical  gyrotron  can  be  used  for 
obtaining  range  resolution  for  isolated  targets  with  roughly  2-3  mm  range 
resolution.  However  one  drawback  is  that  the  target  must  be  stationary, 
and  the  time  to  perform  the  image  is  determined  by  the  time  it  would  take 
the  quasi-optical  gyrotron  to  tune  through  its  entire  tuning  range.  This  is  a 
period  of  several  minutes  as  the  experiment  is  set  up  now.  Several 
approaches  are  possible  to  achieve  much  faster  tuning. 

One  approach  to  rapid  range  resolution  using  the  full  bandwidth  of 
the  device  would  be  a  high  voltage  quasi-optical  gyrotron.  This  would  be 
able  to  tune  through  a  much  large  range  by  changing  the  voltage,  which 
can  be  done  very  quickly.  However  this  experiment  is  only  in  the  concept 
stage  at  this  time.  Another  approach  is  to  use  a  trim  coil  inside  the 
superconducting  coil  of  the  quasi-optical  gyrotron  to  electronically  sweep 
the  magnetic  field.  The  case  of  a  5  kG  trim  field  oscillating  at  0.5  Hz, 
appears  to  be  fairly  straightforward;  the  case  of  a  5  kG  field  oscillating  at  5 
Hz  appears  to  be  possible  but  difficult42.  Half  of  the  period  of  oscillation 
then  defines  the  minimum  time  in  which  the  range  image  can  be  formed. 
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3.  Remote  Sensing  of  Clouds  and  Atmospheric  Particulate  Dispersion 

In  this  section  we  discuss  the  use  of  the  gyrotrons  and  quasi- 
optical  gyrotrons  for  remote  sensing  of  clouds  or  particulates  in  the 
atmosphere.  The  properties  of  many  of  the  different  types  of  clouds 
are  summarized  in  Ref  28. 

A.  Radar  Scatter  from  Clouds 

Since  the  droplet  diameter  in  the  cloud  is  much  less  than  a 
quarter  of  the  wavelength,  the  scattering  is  in  the  Rayleigh  regime 
and  the  cross  section  is  given  by 

o=  7t5|K2|D6A4.  (15) 

where  D  is  the  aerosol  diameter,  X  is  the  wavelength,  and  K  =  (m^- 
1  )/(m^+2)  and  m  is  the  index  of  refraction.  For  water  at  millimeter 
wavelength,  IK^I  is  given  by  about  0.8.  The  strong  increase  of  cross 
section  with  inverse  wavelength  is  a  potential  advantage  of  a 
millimeter  wave  system.  So  is  the  fact  that  a  tightly  focused  beam 
can  be  produced  with  a  relatively  small  antenna.  Since  the  radar 
beam  intersects  the  entire  cloud,  the  radar  cross  section  is  the  cross 
section  of  each  droplet  times  the  number  of  droplets  in  the  beam. 

The  volume  of  the  cloud  examined  is  the  area  of  the  radar  beam 
times  the  range  resolution  distance,  denoted  h.  For  cloud  studies,  the 
radar  reflectance  t]  is  the  usual  parameter  used,  where  "n  = 
IjaiG/4KR2h,  the  cross  section  per  unit  volume,  where  the  summation 
is  over  the  droplets  in  the  volume.  For  a  typical  cloud,  like  those 
studied  by  Lhermitte,  the  droplet  diameter  is  about  lOpm  and  the 
droplet  density  is  about  lOOcm-3,  so  the  reflectivity  is  about  3x10* ^ 2 
cm-1.  With  such  small  reflectances  and  low  power  of  the  transmitter, 
the  observation  of  the  cloud  is  difficult.  For  Lhermitte's  case,  with  a 
range  of  3  km,  a  transmitted  power  of  1  kW,  a  range  resolution 
distance  h  of  6x103  cm,  and  an  effective  aperture  of  4x103  cm^,  the 
received  power  is  about  10*^3  w.  Lhermitte’s  receiver  had  a  noise  of 
about  5x10'^3  so  his  signal  power  is  below  the  noise  power,  and  he 
had  to  integrate  the  signal  over  a  considerable  number  of  pulses  to 
get  good  data.  In  his  case,  he  runs  the  ElO  at  about  10  kHz  pulse 
repetition  rate  and  integrates  for  3  s,  or  about  3x10"^  pulses.  He  finds 
that  the  meteorological  conditions  remain  reasonably  constant  for 
these  3s. 
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Lhermitte  used  an  oscillator,  not  an  amplifier  as  his  microwave 
source.  To  do  the  coherent  signal  integration  and  Doppler  processing, 
he  uses  a  coherent  on  receive  approach  and  corrects  the  phase  of  the 
received  signal  according  to  the  phase  of  the  transmitted  signal. 
Except  for  the  random  phase  from  pulse  to  pulse,  the  EIO  has 
extremely  good  signal  properties.  In  fact,  this  was  the  reason  for 
Lhermitte's  choice  of  an  EIO  over  a  magnetron,  which  could  operate 
at  somewhat  higher  average  power. 


B.  Analysis  of  Radar  Returns  from  Clouds 

Here  we  briefly  examine  the  use  of  a  gyrotron  oscillator  as  a 
source  for  radar  studies  of  clouds  using  the  principles  sketched  out  in 
Section  2.  Thus,  as  in  the  last  section,  we  build  the  data  processing 
scheme  around  the  matched  filter  response,  which  also  optimizes 
signal  to  noise  ratio. 

The  cloud  consists  of  a  number  of  individual  scatterers  denoted 
by  the  index  i.  Each  scatterer  moves  with  velocity  v+5vi,  where  v  is 
the  average  velocity  in  the  range  cell  and  5vi  is  a  random  velocity 
whose  average  value  is  zero.  Each  velocity  is  sufficiently  small  that 
the  scatterers  can  be  regarded  as  stationary  during  the  individual 
radar  pulse.  However  the  velocities  will  be  Doppler  resolved  by 
looking  at  reflections  from  a  long  sequence  of  pulses.  The  return 
signal  from  the  desired  pulse  can  be  constructed  as  in  the  previous 
section.  In  this  case  it  is 

erd(t)expi4)rd(t)  =  li  Si  edt(t-2Ri/c)exp  i<))td(t-2Ri/c) 

X  expi[-(4nifo/c)(Ri +vt+6vit)]  +  Nd(t),  (16) 

where  where  the  summation  is  over  droplets  and  Nd(t)  is  the  noise 
signal  operated  on  by  the  same  operations  that  turn  the  received 
signal  into  the  received  signal  from  the  desired  waveform.  All 
conditions  regarding  the  usable  portions  of  the  spectrum  discussed  in 
the  previous  section  still  apply. 

We  assume  that  the  radar  transmits  P  pulses  spaced  by  a  time 
T,  and  calculate  the  matched  filter  response  for  each  pulse  for  each 
range  cell.  Henceforth  we  will  neglect  the  noise  to  simplify  the 
analysis.  We  assume  that  P  is  large  enough  that  the  signal,  which  is 
proportional  to  P^  is  sufficiently  larger  than  the  noise,  which  is 
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proportional  to  P.  Furthermore,  the  high  power  sources  we  are 
considering  should  be  less  constrained  by  noise  than  was  the  case  for 
Lhermitte.  If  the  desired  pulse  is  Gaussian  with  linear  frequency 
chirp,  the  matched  filter  response  is  also  a  Gaussian  centered  about 
the  center  of  the  range  cell.  We  will  make  the  simplifying 
assumption  that  the  matched  filter  response  is 

Mnp  =  (V7t/2)eo2T£iSiexp-[(47cifo/c)(Ri+(v+5vi)px)]  (17) 

where  the  index  n  denotes  the  range  cell  and  p  the  pulse.  That  is  we 
assumed  each  droplet  is  at  the  center  of  the  range  cell  and  the 
summation  is  over  the  droplets  in  the  range  cell.  Associated  with  the 
index  p  is  the  Fourrier  transform  variable  l%zl?  where  z  takes  on 
integer  values  from  zero  to  P-1.  Then  one  can  define 

Qn(z)  =  SpMnpexp27cipz/P  (18) 

To  obtain  the  reflectance,  one  forms  the  summation  £zlQn(z)l^  =  Qn- 
In  addition  to  the  summation  over  z,  there  are  four  other 
summations  over  i  and  p  and  over  the  analogous  variables  j  and  q  in 
the  complex  conjugate.  The  summation  over  z  is  a  simple  summation 
of  a  geometric  series  and  turns  out  to  be  equal  to  P  if  p=q  and  zero 
otherwise.  We  then  separate  the  summation  over  i  and  j  into  a  sum 
of  two  summations,  one  having  i=j  and  one  having  i;tj.  The  former  is 
a  simple  summation  over  reflectances  over  the  individual  droplets 
and  has  N  terms  in  it  where  N  is  the  number  of  droplets.  The  latter 
is  the  sum  of  0.5N(N-1)  terms,  each  one  having  essentially  a  random 
phase.  The  rms  value  of  a  summation  of  random  phases  is  about 
N,  so  for  each  value  of  p,  the  summation  over  the  i's  and  j's  will  have 
about  the  same  order  of  magnitude.  However  the  summation  over 
the  p  index  greatly  reduces  the  value  of  the  i?tj  summation.  In  this 
summation,  the  sum  over  p  is  also  a  geometric  sum  which  can  be 
done  simply.  The  result  is  that  this  summation  is  very  small  unless 
5vi-5vj<c/4foPx.  Thus  if  the  velocity  spread  of  the  droplets  is 
sufficiently  large,  the  i?*j  sum  is  small  compared  to  the  i=j  sum.  If  fo 
is  94  GHz,  and  Px  =0.1  sec,  the  spread  of  droplet  velocities  only  has 
to  be  large  compared  to  1  cm/s,  a  condition  almost  certainly  satisfied 
for  droplets  in  clouds.  Thus  we  find  that 

Qn=  {0.8P27t7eo4T2d^exp(-2aRn)/16?i6Rn4)XiDi6  (19) 
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where  we  have  used  the  relation  between  s  and  a  for  the  droplets, 
the  expression  for  cross  section  in  Eq.(15),  and  the  fact  that  the 
antenna  gain  is  where  d  is  the  antenna  diameter. 

To  determine  the  average  velocity  of  the  droplets,  one  forms 
the  summation  SzZlQn(z)l^  =  Fn-  An  analogous  calculation  gives  the 
result 


Fn  =  [2PfoVnX/c]xQn  (20) 

where  we  have  explicitly  denoted  the  dependence  of  velocity  on 
range  cell.  Other  moments  of  the  distribution  of  velocities  can  be 
found  by  taking  other  moments  of  Qn.  Equations  (19)  and  (20)  show 
how  the  reflectance  and  relative  velocity  in  a  range  cell  are 
calculated  from  the  backscattered  radar  signal.  Figure  8  is  a 
schematic  of  the  circuitry  for  the  Doppler  cloud  radar. 

C.  High  Power  Radar  Echoes  from  Clouds 

With  this  introduction  to  the  radar  echoes  from  clouds  or 
particulates,  let  us  examine  the  cloud  measurements  Lhermitte  has 
made  and  see  how  they  can  be  enhanced  with  the  use  of  a  high 
power  or  tunable  transmitter.  One  of  the  first  measurements  was 
the  reflection  from  a  thunderstorm  in  the  Miami,  Florida  area.  The 
range  at  which  he  was  able  to  operate  was  about  20-30  km. 

However  since  the  air  at  ground  level  was  quite  humid,  he  was  only 
able  to  observe  the  tops  of  the  clouds.  For  instance  one  of  his 
observations  is  at  an  angle  of  34°  and  a  range  of  20  km,  so  that  he 
observes  an  altitude  of  1 1  km.  He  points  out  that  at  a  30°  elevation, 
the  two  way  absorption  was  12  dB,  whereas  at  a  7°  angle,  it  is  50dB. 
His  transmitter  however  had  an  average  power  of  only  4  W.  A 
gyrotron  transmitter,  operating  with  an  average  power  of  4  kW 
would  be  able  to  examine  the  same  cloud  at  much  lower  elevation. 

One  of  the  easiest  types  of  clouds  to  observe  with  a  94  GHz 
radar  is  high  altitude  cirrus  clouds.  Reference  31  points  out  that 
these  clouds  can  be  observed  with  a  radar  even  when  they  cannot  be 
seen  visually.  The  reason  is  that  a  small  number  of  large  size  ice 
crystals  are  in  the  cloud.  Lhermitte  points  out  that  such  a  cloud  with 
a  single  2  mm  ice  crystal  per  cubic  meter  has  about  35  dB  more 
reflectivity  than  a  fair  weather  cloud  with  the  same  water  content. 
The  tunability  of  the  quasi-optical  gyrotron  could  then  aid  in 
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remotely  sensing  the  distribution  of  ice  crystal  sizes  in  such  clouds. 

As  long  as  the  ice  crystal  diameter  is  less  than  about  a  quarter  of  a 
wavelength,  the  scattering  is  in  the  Rayleigh  regime.  Other  than  the 
dependence  of  the  scattered  power  on  the  wavelength,  there  is  no 
information  concerning  the  distribution  of  sizes.  However  for  larger 
crystals,  there  are  resonances  due  to  Mie  scattering.  Figure  9  shows 
the  backscattering  cross  section  of  water  and  ice  spheres  at  94  Ghz  as 
a  function  of  diameter  (taken  from  Ref.  31).  The  quasi-optical 
gyrotron,  especially  if  it  could  operate  at  both  the  cyclotron 
frequency  and  its  harmonic,  could  irradiate  the  cloud  at  several 
frequencies  between  about  80  and  260  GHz;  the  radiation  could  be 
spread  out  over  several  Mie  resonances,  depending  on  the  size  of  the 
crystals.  The  backscatter  cross  section  as  a  function  of  frequency 
then  gives  information  regarding  the  distribution  of  scatterer  size. 
While  the  size  distribution  cannot  be  exactly  inverted  from  the 
scattered  spectrum,  information  can  be  obtained  as  in  microwave 
propagation  experiments  through  rain'^^ 

The  other  observation  Lhermitte  made  is  of  small  fair  weather 
cumulus  clouds.  A  small  cloud  passes  directly  overhead  at  an 
altitude  of  about  1.5  km;  as  it  goes  by,  the  reflectivity  and  vertical 
velocities  are  measured  by  the  radar  at  a  series  of  range  cells  65  m 
in  depth.  The  cloud  is  overhead  for  about  2  minutes.  Since  each 
vertical  scan  takes  3  s,  there  are  about  40  vertical  scans  in  the  cloud 
and  about  6-8  range  cells.  A  higher  power  transmitter  could 
enhance  these  measurements  in  a  number  of  ways.  First,  in  doing 
the  same  measurement  as  Lhermitte  at  higher  power,  the  radar 
beam  could  be  scanned  from  side  to  side  in  the  cloud  so  that  a  three 
dimensional  image  of  the  cloud  could  be  formed.  Second,  the  range 
cell  could  be  reduced  to  perhaps  10  m  (a  10  MHz  bandwidth)  so  that 
the  cloud  could  be  examined  vertically  in  much  greater  detail.  Third, 
the  higher  power  could  be  used  for  greater  range  to  image  clouds 
that  are  not  directly  overhead.  This  not  only  gives  another  view  of 
the  cloud,  but  it  also  gives  another  component  of  velocity.  For  clouds 
with  a  horizontal  range  of  perhaps  10  km,  the  radar  would  give  a 
Doppler  velocity  in  a  direction  that  is  almost  horizontal. 
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4.  Remote  Sensing  of  Atmospheric  Turbulence 

Here  we  discuss  the  potential  use  of  the  quasi-optical  gyrotron 
as  a  remote  sensor  for  atmospheric  turbulence.  High  power 
gyrotrons  appear  to  provide  a  unique  capability  for  sensing  the  inner 
scale  length  of  the  turbulence  in  an  open  atmosphere  at  fairly  long 
range.  This  is  an  important  property,  because  standard  theories 
relate  it  closely  to  the  strength  of  the  turbulence.  Before  proceeding, 
we  briefly  review  some  properties  of  atmospheric  turbulence.as 
summarized  in  Doviak  and  Zmic37. 

Power  feeding  the  turbulence  is  provided  at  large  scale  length. 

Nonlinear  convection  of  the  fluid  then  generates  shorter  and  shorter 

scale  lengths  until  the  very  short  scale  motions  are  dissipated  by 

viscosity.  This  intermediate  range  of  cascade  is  called  the  inertial 

range.  The  energy  per  unit  mass  of  the  fluid  per  unit  k  is  dE(k)/dk 

and  has  units  m^/s^.  The  power  per  unit  mass  provided  at  small  k  is 

denoted  e,  and  it  has  units  m^/s^.  For  steady  state  turbulence,  e  is 

also  the  power  dissipated  per  unit  mass  at  large  k,  which  is  also 

equal  to  the  power  convected  through  the  inertial  range  in  k  space. 

The  convection  velocity  through  k  space,  dk/dt  ,  is  related  to  E  and  e 

through  the  relation  e  =  (dk/dt)x(dE(k)/dk).  The  assumption  now  is 

that  E(k)  depends  only  on  k  and  e.  Assuming  a  polynomial  relation, 

there  is  only  one  which  is  dimensionally  correct, 

« 

dE(k)/dk  =  A  e2/3k-5/3,  (21) 

where  A  is  a  dimensionless  constant.  Doviak  and  Zmic  use  the  value 
of  e  to  characterize  the  turbulence.  Table  1  gives  their 
characterization.  Generally,  the  spectrum  seems  to  be  borne  out  by 
experiment. 

The  minimum  value  of  k  is  determined  by  the  power  input  at 
large  scale.  The  maximum  value  of  k  is  determined  by  the  viscous 
dissipation.  The  time  for  the  power  in  this  range  of  k  to  dissipate  by 
viscosity  is  given  by  l/vk^,  where  v  is  the  kinematic  viscosity  in 
m2/s.  The  time  the  spectral  energy  remains  in  a  region  of  size  k  is 
given  by  k/(dk/dt).  Equating  these  two  times,  we  find  that  km  ax  is 
proportional  to  £l/4/v3/4.  Doviak  and  Zrnic  give  the  proportionality 
of  1/5,  so 


kmax  =  0.2el/''/v3/''. 


(22) 


At  sea  level,  v  =  10*5  m2/s,  so  for  strong  turbulence,  the  minimum 
distance,  2%lk  is  about  1  cm. 

Scattering  from  clear  air  turbulence  is  caused  by  fluctuations  in 
the  index  of  refraction  which  in  turn  is  caused  mostly  by  fluctuations 
in  water  vapor  content.  To  determine  the  spectrum  of  index  of 
refraction  fluctuations,  the  spectrum  of  a  scalar  quantity  passively 
advected  by  the  turbulent  field  must  be  calculated.37,44.45^  7},e  basic 
assumption  of  all  workers,  however,  is  that  the  passively  advected 
scalar  has  the  same  spectrum  as  the  energy,  so  ln(k)|2  =  Cn^  k*5/3, 
where  Cn^  is  in  units  of  m*2/3.  Doviak  and  Zmik's  characterization 
the  coefficient  Cn^  in  terms  of  the  strength  of  the  turbulence  is 
shown  in  Table  1. 

We  now  consider  the  scattering  from  the  turbulence. 

Assuming  that  5n  is  small,  the  Bom  approximation  for  the  case  that 
the  electric  field  is  evaluated  in  the  far  field  (koD2/r<l,  where  D  is 
the  diameter  of  the  scattering  volume),  gives  the  scattered  electric 
field  at  the  point  r  is  given  by 

Es  =  (ko2Ei/47tr)expikor/d3r’5n(r')exp  i(ki-ks)*r’  (23) 

where  Ej  is  the  incident  electric  field  in  the  scattering  region,  kj  is  the 
incident,  and  ks  is  the  scattered  wave  vector.  The  integral  is  over  the 
scattering  volume,  which  is  defined  as  the  intersection  of  the 
patterns  of  the  transmitting  and  receiving  antennae.  Assuming  that 
the  scattering  volume  is  very  large  compared  to  the  scale  size  of  the 
index  fluctuation,  we  find 

IEs|2=  ko4Ei2VCn2|ki-ksl-^l/3/8r2  (24) 

where  V  is  the  scattering  volume.  In  performing  the  integrals  and 
manipulations,  we  have  made  use  of  the  fact  that  the  spectral 
density  per  unit  k  space  volume  is  given  by  ln(k)|2/47ck2  for  a 
homogeneous,  isotropic  spectrum.  Also,  we  have  made  use  of  the 
fact  that  the  density  of  states  in  k  space  is  V/(27i)3.  If  the  scattering 
angle  is  0,  then  lki-ksl=2kisin0/2. 

For  backscatter  (0=7i)  and  k=1800m-l  (94  GHz),  the  scattering 
wave  vector  is  almost  certainly  larger  than  kmax.  so  there  should  be 
virtually  no  backscatter.  However  there  will  be  scatter  in  other 
directions.  The  use  of  powerful  gyrotrons  at  94  GHz  (and  possibly  at 
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35  GHz  also)  opens  up  the  possibility  of  examining  the  inner  scale 
length  by  oblique  scattering  measurements.  The  inner  scale  length 
could  be  inferred  from  the  angle  at  which  the  scattering  disappears. 
Frequency  tunability,  as  available  with  the  quasi-optical  gyrotron 
could  reduce  the  number  of  receiving  antennae. 

Let  us  consider  a  specific  example  to  show  how  the  availability 
of  94  GHz  gyrotrons  with  over  100  kW  of  power  can  render  possible 
a  scattering  measurement  of  the  inner  scale  length.  The  scattering 
volume  is  then  roughly  LD^  where  L  is  the  length  of  the  intersection 
of  the  two  antenna  patterns  and  we  have  assumed  the  incident  beam 
is  the  more  tightly  focused  of  the  two.  Then  we  relate  the  field 
amplitude  to  the  power  density  and  add  to  Eq.(24)  additional  factors 
accounting  for  the  atmospheric  attenuation.  Then  Eq.(35)  can  be 
rewritten  as 


Pr  =  ko'^AeLPtCn^lkj-ksl'^l/^  expl-a(ri+rs)]/8rs2  (25) 

where  Pr(t)  is  the  incident  (reflected)  power,  Ae  is  the  aperture  of 
the  receiving  antenna,  a  is  the  atmospheric  attenuation,  and  ri  is  the 
range  from  the  incident  antenna  to  the  scattering  region.  Let  us 
choose  a  configuration  in  which  the  scattered  radiation  is  in  the  far 
field.  Although  this  is  not  absolutely  necessary,  it  does  simplify  both 
the  analysis  and  the  interpretation  of  any  data. 

Take  a  transmitting  antenna  of  D=3  m  which  can  focus 
radiation  down  to  a  spot  2  m,  2.5  km  away,  and  a  range  to  the 
receiver  of  10  km.  For  relative  humidity  50%  or  less,  the  attenuation 
is  10  db  or  less.  For  94  GHz  radiation  with  a  3  m  receiving  antenna 
and  an  assumed  inner  scale  length  of  3  cm  (lki-ksl=200),  we  find  that 
Pr/Pi  =  2x10'^^,  assuming  strong  turbulence,  Cn=3xl0’^3^  and  also 
L=D/sin0.  If  the  incident  power  is  10^  W,  the  received  power  is 
about  2x1 0-^2  4  times  the  noise  power  at  300oK.  Integration 

over  a  small  number  of  pulses  could  further  enhance  SNR.  Consider 
the  alternative  of  using  an  EIO  having  a  power  of  1  kW  and  an 
average  power  of  4  W.  Now  the  signal  is  a  factor  of  25  below  the 
noise  level,  and  this  is  for  fairly  strong  turbulence.  It  is  unlikely  that 
the  SNR  could  be  significantly  increased  by  integrating  over  a  large 
number  of  pulses;  at  the  shortest  scale  length  (3  cm),  the  correlation 
time  of  the  turbulence  is  undoubtedly  very  short  (lO'^s  at  3m/s),  so 
the  scattered  pulses  would  not  be  coherent  with  one  another  for  very 
long. 
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5.  Remote  Sensing  of  Humidity 

In  this  section  we  briefly  examine  the  possibility  of  using  the 
high  power  and  tunability  of  the  quasi-optical  gyrotron  to  remotely 
sense  humidity  in  two  ways.  First  we  consider  remote  sensing  on  a 
long  horizontal  path,  second  we  consider  sensing  the  humidity  profile 
of  clouds  by  radar  sensing  in  a  DIAL  mode. 

A.  Horizontal  Path  Average  Humidity  Measurement 

As  pointed  out  in  Ref.  28,  many  attenuation 
measurements  have  been  made,  both  in  the  laboratory  and  in  the 
field.  These  are  difficult  measurements  due  to  the  necessity  of 
keeping  the  electronics  sufficiently  stable.  Generally  the 
experiments  use  cw  sources  chirp  the  frequency.  The  measurements 
we  discuss  here  are  basically  single  or  multi  .frequency 
measurements  with  pulses  several  microseconds  long.  The  high 
power  and  short  pulse  length  are  advantageous  because  the 
conditions  of  the  atmosphere  will  be  constant  for  one  series  of 
measurements.  The  quasi-optical  gyrotron  provides  a  unique 
capability  in  that  it  can  be  tuned  (with  second  harmonic  operation)  to 
the  center  of  the  absorption  peak  at  180  GHz,  and  it  has  enough 
power  to  propagate  long  distances  in  the  atmosphere. 

Operating  the  quasi-optical  gyrotron  at  several  modes  near  180 
GHz  allows  several  frequencies  within  the  absorption  peak  to  be 
examined.  Furthermore,  not  only  absorption,  but  also  phase  shifts, 
which  may  be  easier  to  calibrate,  can  be  measured.  By  comparing 
with  dispersive  delay  times  Xd,  also  shown  in  Fig.  6,  water 
concentration  can  be  inferred.  The  phase  shift  A<P  is  given  in  terms 
of  the  dispersive  delay  time  by  AO  =  27tfxd.  If  the  receiver  is  at 
range  R  from  the  transmitter,  and  the  transmitted  pulse  at  frequency 
fi  (one  of  the  modes  of  the  oscillator)  is  given  by  Ei(t)cos[2jifit-i-(|)i(t)], 
the  received  pulse  will  be  given  by 

Er  =  Eri(t)cos  [27tfit-t-<t)ri(t)]  =  Ei(t-R/c)exp(-aiRcos)[2jifit+())ri(t)],  (26) 

where  ai  is  the  attenuation  and  <;>ri(t)  =27tfiR/c  +  <l)i(t-R/c)  +  AOj  The 
E's  and  <t>’s  can  be  measured  from  the  transmitted  and  received 
signal,  as  discussed  in  Section  2.  Assuming  that  the  range  is  not 
known  accurately  to  a  fraction  of  a  wavelength,  it  can  be  eliminated 
between  the  various  modes.  Specifically,  we  find 

fj[<t)ri(t)-(t)i(t-R/c)]  -  fi[(|)rj(t)-<!)j(t-R/c)]  =  fjAOi-fiAOj  (27) 
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so  all  phase  shifts  can  be  found  in  terms  of  a  single  reference  shift. 
Thus  from  measurements  of  the  phase  and  attenuation,  the  relative 
humidity  can  be  inferred. 

B.  The  Humidity  Profile  of  Clouds. 

Another  application  is  the  remote  sensing  the  humidity  profile 
of  clouds  with  the  use  of  radar  backscatter  at  two  different 
frequencies.  The  droplets  serve  as  the  scatterers  and  the  relative 
attenuation  between  the  different  range  cells  in  the  cloud  can  give 
information  to  the  humidity  profile  in  the  cloud.  As  is  apparent  from 
Eq.(19),  the  reflectivity  from  a  range  cell  in  a  cloud  depends  on  the 
number  of  scatterers  as  well  as  the  attenuation  between  the 
transmitter  and  the  reflecting  region  of  the  cloud.  Consider  the 
return  from  the  range  cell  in  the  cloud  which  is  nearest  the 
transmitter.  At  a  single  frequency,  only  the  product  of  the  scattering 
and  the  attenuation  from  the  transmitter  to  the  bottom  of  the  cloud 
can  be  inferred.  However  by  using  two  frequencies,  the  attenuation 
and  the  scattering  can  each  be  obtained  as  long  as  the  attenuation  is 
known  as  a  function  of  frequency  (as  given  in  for  instance  Fig.6).  In 
doing  so,  one  must  also  account  for  the  fact  that  the  summation  over 
droplets  may  have  a  frequency  dependence  due  to  the  fact  that  the 
geometry  of  the  scattering  volume  may  also  depend  on  frequency  for 
a  specified  antenna.  By  stepping  up  in  the  cloud  from  one  range  cell 
to  another,  the  humidity  profile  can  be  obtained  within  the  cloud. 

The  quasi  optical  gyrotron,  which  has  the  capability  of  both  high 
power  and  tunability,  could  provide  a  unique  source  source  for  such 
a  measurement  by  examining  the  cloud  at  first  90  and  then  130  GHz. 
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6.  Remote  Sensing  of  Trace  Impurities  in  the  Upper  Atmosphere 

We  now  look  into  whether  the  tunability  of  the  quasi-optical 
gyrotron  can  be  used  to  detect  trace  elements  other  than  water.  At 
sea  level,  this  appears  to  be  an  extremely  difficult  measurement  due 
to  the  fact  that  the  pressure  broadening  of  the  absorption  lines  is 
substantial  (typically  about  2-3  GHz),  and  the  impurity  concentration 
is  very  low.  Distinguishing  the  impurity  absorption  from  the  water 
vapor  absorption,  which  dominates  it  by  many  orders  of  magnitude, 
does  not  appear  to  be  a  simple  task.  However  for  impurities  in  the 
upper  atmosphere,  the  pressure  broadening  is  much  less,  as  is  the 
amount  of  water  vapor  to  interfere  with  the  measurement. 

The  idea  is  then  to  bounce  the  radiation  at  a  frequency  near 
the  absorption  line  off  an  orbiting  satellite,  detect  the  backscattered 
radiation,  and  examine  the  absorption  and  phase  shift  as  a  function 
of  frequency  to  determine  the  concentration  of  the  impurity.  (One 
might  think  that  at  high  altitude,  the  molecule  would  spontaneously 
re-emit  before  it  is  collisionally  deexcited,  so  that  a  satellite  would 
not  be  necessary.  However  the  time  for  spontaneous  reemission  for 
rotationally  excited  molecules  varies  from  about  10“^  to  10^  s,  so  this 
is  not  the  case^^.)  Measurements  like  this  are  currently  performed 
by  radiometry,  whereby  the  thermal  spectrum  of  the  atmosphere  is 
measured  and  related  to  the  trace  elements.  The  scattering  scheme 
we  discuss  here  would  appear  to  have  several  potential  advantages 
over  radiometry.  First  of  all,  the  effect  is  doubled  because  two  way 
propagation  doubles  the  path  length.  Secondly,  phase  measurements 
as  well  as  amplitude  measurements  now  become  possible.  Finally 
because  the  power  of  the  millimeter  wave  sources  is  large,  the 
measurements  give  snapshots  of  the  atmospheric  concentration  in 
very  short  times,  milliseconds  to  seconds  depending  on  the 
concentration  and  antenna  size.  This  contrasts  with  the  many  hours 
or  even  days  it  takes  to  do  an  analogous  measurement  by  ground 
based  radiometry. 

A.  Trace  Element  Determination  by  Radiometry 

In  radiometry,  the  thermal  spectrum  of  millimeter  wave 
radiation  is  detected  and  related  to  the  concentration  of  the 
element. These  measurements  can  be  absorption  measurements, 
for  instance  measuring  the  absorption  of  sunlight  at  frequencies  near 
the  absorption  line.  Also  the  thermal  emission  can  be  studied,  and 
this  is  what  we  concentrate  on  here.  For  the  atmosphere  which 


radiates  as  a  black  body  (if  it  were  optically  thick)  the  radiation 
temperature  is  related  to  the  intensity  of  radiation  by  the  Rayleigh- 
Jeans,  or  actually  the  Planck  law.  If  the  radiation  is  emitted  from  a 
region  with  temperature  T  and  s  is  the  path  length,  then  the 
radiation  temperature  Tr  at  the  receiving  antenna  is  given  by 

Tr  =  T(s)a(s)exp-t(s),  t(s)  =  lha(s’)ds’  (28) 

where  a(s)  is  the  sum  of  the  absorption  coefficient  of  all  gases  to  be 
considered  and  x  is  the  opacity.  For  a  propagation  medium  which  is 
optically  thin,  as  is  the  case  for  millimeter  waves  in  the  upper 
atmosphere,  the  presence  of  the  gas  of  interest  gives  rise  to  a  slight 
increase  in  radiation  temperature  at  the  frequencies  near  the 
absorption  line. 

The  upper  atmosphere  i  typically  at  a  temperature  of  between 
about  200-300OK.  Thus  for  an  optically  thin  propagation  medium, 
the  ratio  of  antenna  temperatu:c  due  to  ‘he  impurity  to  the  source 
temperature  (that  is  the  upper  atmosphere)  is  roughly  the  opacity. 
For  ozone  at  110.8  GHz,  in  a  standard  atmosphere,  the  zenith  opacity 
is  about  6%  and  the  antenna  temperature  is  about  10°  K  at  the  line 
center. if  the  line  shape  is  Lorentzian,  as  is  characteristic  of  a 
pressure  broadened  line,  the  temperature  a^  a  function  of  frequency 
reflects  this  line  shape.  The  line  width  is  a  function  of  altitude.  If 
sufficient  resolution  of  temperature  as  a  function  of  frequency  is 
available,  the  profile  of  ozone  density  as  a  function  of  altitude  can 
often  be  inverted.4  8 

This  temperature  due  to  the  impurity  is  generally  small 
compared  to  the  atmosphere  temperature.  In  radiometry,  the 
relative  temperature  fluctuation  through  a  channel  of  bandwidth  B  is 
given  by  (BY)- 1/2  where  Y  is  the  integration  time,  so  small 
temperature  changes  due  to  an  impurity  can  be  resolved  by  long 
integration  time.  Parrish  et  al49  have  recently  studied  trace  element 
concentrations  in  the  upper  atmosphere.  These  are  all  related  to  the 
ozone  cycle,  and  as  such,  the  altitudes  are  about  40  km,  so  the 
pressure  broadening  is  about  10-20  MHz.  They  found  that  at  277 
GHz,  ozone  had  a  temperature  of  8°K  and  the  integration  time  was  an 
hour  to  achieve  SNR  of  600.  They  also  found  that  CIO  had  a 
temperature  about  0.06°  at  279  GHz.  To  achive  SNR  of  25  took  36 
hours  over  9  days.  The  weakest  line  they  were  able  to  detect  was 
the  266  GHz  line  of  HO2.  This  line  had  a  temperature  of  about 
0.01 5°K  and  took  55  hours  to  resolve  the  line. 
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B.  Ozone  Detection  by  Backscatter  from  an  Existing  Satellite 

We  now  consider  the  use  of  the  quasi-optical  gyrotron  as  a 
tunable  source  for  active  detection  of  ozone.  We  consider  the  110.8 
GHz  absorption  line;  a  frequency  at  which  the  quasi-optical  gyrotron 
can  easily  be  tuned.  To  provide  a  two  way  path,  we  exploit  reflection 
from  an  existing  satellite.  There  are  many  satellites  and  pieces  of 
space  debris  whose  low  earth  orbits  are  known  to  high  precision. 

As  we  will  see  however,  not  all  are  suitable  as  reflectors;  to  be  a 
viable  reflector,  the  size  of  the  satellite  must  be  quite  small. 

In  active  probing  of  the  upper  atmosphere,  both  amplitude  and 
phase  are  available  to  analyze.  Here  we  consider  only  the  phase, 
since  the  analysis  of  the  amplitude  is  no  different  in  principle  from 
radiometry.  For  a  Lorentzian  line  shape  centered  at  (Oq,  the  total 
complex  attenuation  and  phase  shift,  H  for  the  two  way  propagation 
is  given  by 

E  =  Tv((v  +  i(a)-o)o))/(v2  +  (©-cdo)^)},  (29) 

where  t  is  the  two  way  opacity  at  the  line  center,  and  v  is  the 
pressure  broadening,.  At  the  altitude  in  question,  vHti  is  about  10 
MHz.  Thus,  for  ozone,  whose  two  way  zenith  opacity  is  12%,  and 
whose  opacity  at  45°  is  about  20%,  the  total  phase  shift  as  one 
crosses  the  line  could  be  as  high  as  100-15°.  Analogous  to  Eq.  (26),  if 
the  incident  wave  is  at  frequency  fj  the  phase  of  the  returned  signal 
is  given  by 

<!)ri=  27rR(2fi+5fdi)/c  +  <|)ti(t-2R/c)  +  A<I)i(fi) 

+  AOi(fi+5fdi)  +0  +5<J>(fi,p)  (30) 

For  a  satellite  of  known  orbit,  the  Dopp’er  shift  can  be 
calculated  very  accurately,  and  can  be  as  high  as  a  few  MHz  for  100 
GHz  incident  radiation.  In  Eq.(30),  A<l>i(fi)  is  the  ozone  generated 
phase  shift  at  frequency  fi  on  the  way  up,  AOi(fi+8fdi)  is  the  ozone 
generated  phase  shift  of  the  Doppler  shifted  wave  on  the  way  down, 
0  is  the  phase  shift  from  all  other  components,  determined  from  for 
instance  a  standard  atmosphere  model,  and  5<I)(fi,p)  is  the  phase  shift 
generated  by  the  satellite,  where  p  denotes  pulse  number.  The 
satellite  generated  phase  shift  varies  with  frequency  because 
different  frequencies  have  different  scattering  centers;  it  varies  with 


p  because  the  satellite  might  change  its  orientation,  due  to  rotation, 
during  the  pulse  train.  The  difficulty  is  that  5d>(fi,p)  is  basically 
unknown.  One  can  only  perform  the  measurement  if  5<l)(fi,p)  is 
independent  of  both  fi  and  p. 

The  satellite  will  be  considered  to  be  a  set  of  discrete 
scatterers,  with  index  m,  centered  at  R  and  having  position  5Rm  from 
the  center.  It  is  the  5Rm’s  which  contribute  the  essentially  random 
nature  of  the  phase.  If  the  frequency  difference  between  the  waves 
in  the  spectrum  is  denoted  5f,  then  the  contribution  to  the  phase 
shift  is  4ji6R8f/c,  where  for  convenience,  we  have  deleted  subscripts. 
As  long  as  this  is  small  compared  to  the  phase  shift  we  are 
attempting  to  measure,  it  can  be  neglected.  This  then  puts  a  limit  on 
both  the  size  of  the  satellite  and  also  on  the  frequency  spread.  For 
ozone  at  40  km  altitude,  assume  a  width  of  20  MHz..  At  a  45° 
inclination  angle,  the  phase  shift  we  are  measuring  is  about  0.2 
radians.  This  means  that  the  satellite  or  debris  radius  must  be  small 
compared  to  about  30  cm.  Also,  the  phase  might  change  pulse  to 
pulse  because  of  the  rotation  of  the  satellite,  so  the  total 
measurement  time  must  be  small  compared  to  a  rotation  period  of 
the  satellite,  typically  tens  of  seconds.  Thus,  a  small  existing  satellite 
or  piece  of  space  debris,  which  is  not  rotating  violently,  could  serve 
as  a  reflector  for  rapid,  ground  based  ozone  measurement. 

We  now  consider  what  it  would  take  to  do  an  ozone 
measurement  at  a  transportable  facility.  If  we  consider  the  antenna 
to  be  3  m  (a  gain  of  70  dB),  a  bandwidth  which  is  pulse  time  limited 
(10^  s'l),  a  transmitted  power  of  1  MW,  a  range  of  300  km,  a  5  cm 
radius  spherical  scatterer,  and  a  receiver  temperature  of  lOOQoK,  we 
find  that  SNR  of  10^.  This  corresponds  to  a  phase  resolution  in  a 
single  pulse  of  about  (SNR)' 1/2.  if  20  pulses  within  a  time  of  a  tenth 
of  a  second  to  a  second  are  transmitted  within  the  20  MHz 
interrogation  bandwidth,  the  phase  shift  within  the  absorption  line 
can  be  resolved  and  the  ozone  density  can  be  inferred. 

C.  Measurement  of  Trace  Impurity  Density  with  a  Specially  Designed 

Satellite 

The  satellite  generated  phase  shift  can  be  greatly  reduced  or 
eliminated  by  the  use  of  a  specially  designed  satellite.  Specifically, 
we  consider  the  use  of  a  spherical  satellite,  whose  cross  section  is 
independent  of  orientation,  and  whose  cross  section  as  a  function  of 
frequency  can  be  accurately  calculated  and  measured.  For  trace 


elements  such  as  CIO  or  HO2,  discussed  in  Ref.  49  ,  the  total  phase 
shift  across  the  absorption  line  is  10‘3-1 0*4  radians.  Thus  the 
satellite  must  be  smooth  to  about  this  fraction  of  a  wavelength. 
However  since  we  are  considering  millimeter  waves,  this  fraction  of  a 
wavelength  means  an  optical  quality  reflector,  a  standard  capability. 
We  consider  specifically  a  20  cm  radius  spherical  reflector.  In  one 
rocket  launch,  a  number  of  these,  perhaps  ten  or  so,  could  be  placed 
into  orbit  to  provide  nearly  continual,  permanent  coverage. 

We  now  consider  reflection  from  one  of  these  specially 
designed  satellites  to  measure  trace  impurity  density.  For  instance 
the  242  GHz  CIO  absorption^O  line  would  be  observable  with  the 
quasi-optical  gyrotron  at  the  second  harmonic.  It  has  a  limb 
temperature  of  few  degrees  Kelvin,  and  a  ground  temperature  of 
about  O.O50K  looking  near  the  horizontal.  Since  an  active 
measurement  has  double  the  path  length,  we  assume  that  at  an  angle 
of  450  to  the  zenith,  a  value  of  about  O.O30K,  or  a  phase  shift  of  about 
3x10*4  radians.  The  total  two  way  absorption,  from  sea  level,  at  240 
GHz,  at  an  inclination  of  45°  to  the  zenith  is  less  than  10  dB  for 
relative  humidity  less  than  50%.  If  a  high  altitude  is  chosen  for  the 
measurement,  the  absorption  is  less  still.  We  consider  first  the  use  of 
a  satellite  tracking  station,  with  an  antenna  diameter  of  20  m. 
Consider  a  receiver  temperature  now  of  lOO^K,  a  total  atmospheric 
attenuation  of  10  dB,  and  a  radiated  power,  assumed  to  be  100  kW  at 
the  second  harmonic.  Then  for  a  10  ps  pulse,  SNR  is  about  2x10^  . 

Ten  pulses  will  then  give  sufficient  SNR  to  accurately  measure  a 
phase  shift  3x10*4  radians.  This  time  is  so  short  that  it  could  even  be 
possible  to  do  the  measurement  over  a  longer  time  with  a  3m  dish  in 
a  portable  facility.  However  the  satellite  would  have  to  be  accurately 
tracked  over  the  measurement  time.  To  conclude,  the  tunability  of 
the  quasi-optical  gyrotron  gives  a  capability  to  rapidly  measure 
upper  atmosphere  ozone  concentrations  with  existing  equipment. 

With  specially  designed  satellites  and  accurate  satellite  tracking,  it 
could  give  rapid  measurement  of  trace  impurities  at  much  lower 
concentration. 
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Table  1  —  The  Characterization  of  Turbulent  Strength 


Turbulence 

Cn2  in-2/3 

e  m2/s-5 

Weak 

6x10-17 

3x10-3 

Moderate 

2x10-15 

8.5x10-3 

Severe 

3x10-13 

6.8x10-2 
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Fig.  6  —  Attenuation  (dB/km)  and  phase  delay  (ps/km)  for  atmospheric 
propagation  as  a  function  of  frequency  and  1  lidity 
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Circuit  model  for  the  qiiasi-optical  gyrotron  fur  super  range  resolution.  After  the  A/D  converter, 
the  circuit  elements  are  computer  programs  and  the  data  are  analyzed  off  line. 
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